Connection between supernova shocks, flavor transformation, and the neutrino signal 
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We show that the supernova neutrino signal in terrestrial detectors could give characteristic sig- 
natures as the shock propagates through the regions where matter-enhanced neutrino flavor trans- 
formation takes place. This effect could allow existing and proposed supernova neutrino detectors 
to provide new insights into the problem of shock re-generation and propagation. Alternatively, 
such observations could allow unique insights into the neutrino mass and mixing matrix. 



PACS numbers: 14.60.Pq, 12.15.Ff, 97.10.Cv, 97.60.Bw 

In this letter we point out a link between the expected 
neutrino signal from a core collapse supernova event on 
the one hand, and neutrino mass/mixing physics and su- 
pernova shock re-heating and propagation physics on the 
other. The crux of the supernova problem is understand- 
ing how the initial core bounce shock (which stalls at 
time post core bounce i p b ~ 0.1s) becomes re-energized 
by neutrino heating (at i p b < Is) . Conventional ob- 
servations (optical, radio, and x-ray) of supernova events 
and remnants give little direct information on the re- 
heating history and early propagation of the shock. Here 
we demonstrate that the neutrino signal may show char- 
acteristic time structure in energy /luminosity stemming 
from shock propagation across the (resonance) regions 
where matter-enhanced neutrino flavor transformation 
takes place. In principle, observations of such feaures 
could provide insight into the shock strength and loca- 
tion as functions of time. The experimentally-determined 
mass and mixing properties for the active neutrinos sug- 
gest that these flavor transformation-generated features 
coincidentally occur when the neutrino luminosity is 
high. Detection of such shock-related time structure in 
the supernova neutrino signal, in turn, could yield infor- 
mation on neutrino mass and mixing parameters compli- 
mentary to that derived from experiment. 

Solar and atmospheric neutrino experiments have re- 
vealed that neutrinos have vacuum masses and that they 
may mix with large (perhaps near maximal) mixing an- 
gles H §). The LSND experiment j§ can be interpreted 
as evidence for additional vacuum neutrino mixing in the 
f± v e channel with a characteristic mass-squared dif- 
ference Sm 2 ~ 1 e V 2 . This result will be checked in 
the near future by the MiniBooNE experiment ||. If 
this signal holds up, then in combination with the so- 
lar and atmospheric neutrino anomalies we are forced 
to adopt new physics: either the introduction of a light 
singlet ("sterile") neutrino which mixes with the active 
species; or CPT-violating schemes in which neutrinos and 
antincutrinos have different masses and mixing proper- 
ties |1. If the former, then there are two mass/mixing 
arrangements which have been advanced as giving si- 
multaneous solutions for all the neutrino anomalies: the 



"3 + 1" scheme wherein the fourth (heaviest) neutrino 
mass state could be mostly singlet in character; and the 
"two-doublet" scheme in which a lower mass doublet is 
involved in the solution of the solar neutrino problem and 
the higher mass doublet is invoked for the atmospheric 
anomaly and the LSND signal. In any case, our neutrino 
flavor transformation/shock propagation features in the 
supernova signal will arise for the pure active (3 neutrino) 
scenario as well as for those involving sterile neutrinos. 

Matter enhancement of neutrino flavor conversion in 
the supernova will be engendered by these neutrino 
mass schemes via the Mikheyev-Smirnov-Wolfenstein 
(MSW) mechanism § §. MSW resonances in the region 
above the hot proto-neutron star could alter the energy 
spectra and fluxes of the expected standard supernova 
neutrino signal. Neutrino flavor conversion effects on the 
supernova neutrino signal have been considered [jl(| H 
previously, but without consideration of shock effects. 

In fact, the shock represents a major perturbation on 
the pre-supernova density and electron fraction (Y e ) pro- 
files. Furthermore, the shock is expected to be moving 
through the region above the neutron star where neutri- 
nos with typical energies (E v ~ 10 MeV) will be resonant 
at the times (i p b ~ 0.1s to ss 15 s) when neutrino fluxes 
are large enough to be potentially detectable. Neutrino 
flavor conversion efficiency at resonance depends on the 
density/electron fraction gradient. Since the shock front 
is very steep, neutrino flavor conversion can be "shut off" 
as the shock passes through an MSW resonance region. 
Subsequently, as the shock moves beyond the resonance 
region, flavor conversion could be reestablished, albeit at 
a smaller radius, likely corresponding to a higher den- 
sity gradient. The net result could be that, for example, 
V{i,t v<s flavor conversion is "interrupted" for a time of 
order the shock (front and rarefaction zone) propagation 
across resonance regions. This can give dips or steps in 
neutrino energy/count rate in terrestrial detectors. 

Core collapse supernova simulations including numeri- 
cal treatments of hydrodynamics and neutrino transport 
are not always in agreement. However, there are essen- 
tial features of shock physics that are common to these 
models: a rapid (instantaneous) density jump, followed 
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FIG. 1: Adopted density profile for progenitor and shock wave 
as a function of radius at various times post-bounce t p h- 



by a rarefaction region where the density can drop well 
below that at the leading edge of the shock. We have 
drawn on results from several calculations |fl2| to produce 
a generic time-dependent fit to the density profile imme- 
diately above the neutron star and through the region of 
the shock. The density profiles at several times are shown 
in Fig. 1. We have employed ap~ 1/r 24 profile for 
both the progenitor star and the "neutrino-driven wind" 
which forms above the neutron star at i p b > 2 s. (The re- 
sults are similar if we use a p ~ 1 /r 3 profile as suggested 
by the progenitor model of Ref. Q.) We have steepened 
the shock front in our density profiles relative to that 
suggested by numerical calculations. Shock fronts in hy- 
drodynamics calculations may be softened by numerical 
techniques, leading to unphysically large neutrino flavor 
conversion probabilities there. 

We take the neutrino energy distribution func- 
tions at the neutrino sphere to be Fermi-Dirac black 
bodies with small chemical potentials, f{E v ) — 
(l/T^F 2 (ri Ua ))£?/(exp (E v /T Ua - r) v<x ) + l), for neu- 
trino species v a with a = e, /a, t (assuming no singlet neu- 
trino flux from the neutron star surface). Here average 
neutrino energies are given in terms of neutrino sphere 
temperatures by (E Va ) = T,y a F 3 {r}v a ) / F 2 {'q Va ) with F k 
the standard relativistic Fermi integral of order k. For 
illustrative purposes, we have adopted neutrino degener- 
acy parameters rj Va — 3 for all neutrino species and have 
taken (E v ^) = (% (T ) = 26.9 MeV, (E Pe ) = 16.0 MeV, 
and (E Ue ) = 10.0 MeV for all times Q. In fact, neutrino 
degeneracy parameters and average energies will evolve 
with time. At early times it may be that all neutrino 
species have similar energy spectra, but at late times 
when the core has substantially de-leptonized, they may 
be different as above. 

The energy differential fluxes for each neutrino species 
are given in terms of the luminosities L Va , neutrino 
sphere radius i?„, and radius from the center of 



the neutron star r by dcj){E v ) w (L„ Q /27ri? 2 )[(l — 
y/1 - Rl/r 2 )/{E Ua )]f{E v )dE v , We take°the liberated 
gravitational binding energy as 3 x 10 53 ergs; and neu- 
trino luminosities to be the same for all active species 
with an exponential decay time constant of t = 3 s. 

Neutrino flavor amplitude evolution here is calculated 
as in Rcf.s |], in an assumed succession of two-by 
two transformations: 
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where a a and a/3 are the amplitudes for the neutrino to 
be flavor v a or i/p, respectively; A = Sm 2 /2E V with 
Sm 2 the appropriate neutrino mass-squared difference. 
The neutrino potentials which drive flavor evolution are: 
A = ±v / 2G , F^t>}'e, where rib is the baryon number den- 
sity; B and B s are the flavor diagonal and off-diagonal 
neutrino-neutrino forward scattering potentials, respec- 
tively. In the calculations presented here we neglect the 
neutrino-neutrino forward scattering contributions to the 
potentials. While justified where the neutrino fluxes are 
low; at early times, near the neutron star, or where the 
entropy is high, these terms (as well as feedback on Y e ) 
can be important in determining MSW flavor transforma- 
tion probabilities ]l5|) . These effects can alter our time- 
dependent spectral features in some cases 

We discuss the supernova neutrino signal resulting 
from four distinct neutrino mass and mixing schemes: 
(1) 3-active neutrinos only, with a "normal" mass hier- 
archy; (2) "3 + 1" with m V4 ^ m V3 ; (3) the two doublet 
scheme; and (4) a CPT-violating scheme with the neu- 
trinos as in (1) but with an inverted mass hierarchy in 
the anti-neutrino sector producing additional resonances. 
We find that shock effects are important in setting the 
neutrino signal in all of these neutrino mass schemes. 

For the 3-active neutrino scheme and typical density 
profiles across and behind the shock front, there will arise 
a sequence of MSW resonances (each in the f M)T ^ v e 
channel): up to three resonances associated with the 
atmospheric mass-squared difference scale (taken to be 
8m 2 = 3 x 10~ 3 eV 2 ); and up to three associated with the 
solar 5m 2 -scale. (For all cases we take the Large Mixing 
Angle solar solution and adopt Sm 2 — 3.7 x 10~ 5 eV 2 
and sin 2 28 & = 0.8.) The resonance density is (pY e ) res f» 
6.55 x 10 6 gcm -3 (5m 2 cos 26 /E v where Sm 2 is in eV 2 and 
E v is in MeV. For demonstration purposes, we take 
#13 = 0.07. It should be noted that several features in 
the neutrino signals are sensitive to the choice of #13, 
possibly allowing a probe of this parameter. 

The "3 + 1" scheme has an additional set of higher- 
density resonances associated with the LSND 8m 2 -scale 
in the v e ^ v s channel. In this case, Y e is replaced 
by 3(Y" e ^ l/3)/2. (For illustrative purposes we adopt 
5m 2 = 1.0 eV 2 and sin 2 26> LSND = 3.5 x 10~ 3 for all 
LSND-scale two-by-two resonances, where 6>lsnd should 
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be understood to be an effective two-by-two vacuum 
mixing derived from the appropriate 4-neutrino unitary 
transformation.) For the two-doublet scheme we follow 
Ref.s [ jl7[ and exploit the (near) maximal f± v t mix- 
ing to treat flavor evolution above the neutron star as 
two sequential two-by-two level crossings: i^ )T f± v a fol- 
lowed at larger radius (lower density) by v e ^± v ix,t- At 
even larger radius there can be a matter-enhanced tran- 
sition at the solar 8m 2 scale. Given the shock profiles as 
in Fig. 1, each level crossing could give up to three res- 
onance regions. We take an inverted mass hierarchy for 
the antineutrino sector of case (4) . The LSND <5m 2 -scale 
is associated with the v e ==± P MiT resonance in this case. 

We calculate the evolution of the neutrino flavor states 
through our adopted density profile (including the time 
dependent position of the shock) by numerically integrat- 
ing Eq. [l] using an adaptive Runge-Kutta procedure or 
the semi-classical Landau-Zener approach |1| (with sub- 
sequent projection onto vacuum flavor eigenstates) where 
applicable. As the shock moves out and overtakes the 
MSW resonance regions appropriate for each neutrino 
mass scheme, the neutrino energy spectra for v^^ and v e 
will evolve as a result of shock modification of the den- 
sity profile. In particular, the rarefaction zone behind 
the shock may introduce additional level crossings (reso- 
nances) that would be absent in the smooth, monotonic 
density run of the progenitor star. 

Fig. 2 shows charged current (CC) event rates for v e 
and v e and total neutral current (NC) rates on deu- 
terium |l9| as functions of time from a 10 kpc distant 
supernova for the four mass/mixing schemes discussed 
above. These rates are calculated for 0.77 kT of D2O, 
similar to the fiducial mass in the SNO detector. Ratios 
of CC to NC rates are plotted in Fig. 3. The variation 
in time (and associated model-dependent uncertainties) 
of the neutrino emission rates from the core are thereby 
suppressed. Statistical fluctuations in the signals are ne- 
glected, and we take Y e — 0.5 in this simplified example. 
All four mass/mixing schemes show deviations in time 
from a simple exponential for the event rates, or from a 
flat line for the ratios. The two doublet scheme shows a 
dip in the v e CC event rate centered at approximately 1 s. 
When integrated over the full width, this dip represents 
a 33% reduction out of an estimated 58 events over the 
same time period without this effect. The CPT violating 
scheme exhibits a 37% reduction out of ^46 events in 
the v e CC rate for a dip also centered at about 1 s. The 
3-active, "3 + 1," and CPT violating cases show similar 
dips in the v e CC rate in the range 5 s < i p b < 11 s 
representing about a 20% reduction out of ^60 events. 

While the timing of the structures in the CC signals 
are dependent upon the assumed shock model, we ex- 
pect that such structures will be a generic feature of the 
neutrino burst from a core-collapse supernova. NC rates 
show time structure for the cases which include a ster- 
ile neutrino, e.g. two doublet and "3 + 1", though at 
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FIG. 2: Event rates for v e CC (solid), u e CC (dotted), and 
NC (dot-dash) in SNO as functions of time for: (a) 3-active, 
(b) "3 + 1," (c) two doublet, and (d) CPT violating schemes. 
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FIG. 3: Ratios of u e (solid) and v„ (dotted) CC to NC event 
rates for SNO. Neutrino mixing schemes are (a) 3-active, (b) 
"3 + 1," (c) two doublet, and (d) CPT violating. 



low amplitudes for the parameters discussed here. The 
3-active and CPT-violating cases would show no shock- 
related time dependence in the NC signal since the NC 
interactions are "flavor blind." Improved signal statistics 
could be obtained by considering detectors with substan- 
tially increased target mass. The use of light/heavy water 
and -CH2 as supernova neutrino detection media are dis- 
cussed in Ref.s |^0|, ^IJ. The proposed use of high atomic 
number materials such as lead llll W\ could significantly 
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imp rove the detectabilitv of the structures in the v„ CC 
sign kls. 
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For the results shown in Figs. 2 and 3, we have as- 
sumed all resonances are adiabatic except at the shock 
front, including those resonances which appear momen- 
tarily straddling the rarefied region behind the shock. 
Even in the cases where the two resonance regions in the 
rarefaction zone overlap, our numerical treatment shows 
that the flavor transformation from the first (smaller ra- 
dius) resonance is effectively nullified by the second. If 
one assumes that the density profile behind the shock 
includes significant stochastic fluctuations, then this pre- 
scription is changed somewhat. In the case where "flavor 
depolarization" is operative |2^| , the resonances behind 
the shock would yield an equal mix of flavor states. Such 
a treatment yields CC signals similar to that above, but 
with the late count rate reduced. For the situation where 
the density fluctuations or gradients behind the shock are 
sufficiently strong to preclude adiabatic transformations 
entirely, the late time CC ratios are reduced even fur- 
ther, and the "dips" become steps. (The signal from 
SN1987a has been analyzed with respect to the CPT- 
violating, inverted £/-mass scheme |24|] without including 
shock effects, with a "tension" noted between the en- 
ergy/count rate for the detected events and that expected 
from SN models. If flavor depolarization or non- adiabatic 
level crossings occur behind the shock, this apparent dis- 
crepency may be reduced by the effects described here.) 
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